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ABSTRACT 

We present a statistical analysis of the X-ray bursts observed from the 2002 
June 18 outburst of the Anomalous X-ray Pulsar (AXP) IE 2259-1-586, observed 
with the Proportional Counter Array (PCA) aboard the Rossi X-ray Timing 
Explorer. We show that the properties of these bursts are similar to those of 
Soft Gamma- Repeaters (SGRs). The similarities we find are the burst durations 
follow a log-normal distribution which peaks at 99 ms, the differential burst 
fluence distribution is well described by a power law of index —1.7, the burst 
fluences are positively correlated with the burst durations, the distribution of 
waiting times is well described by a log-normal distribution of mean 47 s, and the 
bursts are generally asymmetric with faster rise than fall times. However, we find 
several quantitative differences between the AXP and SGR bursts. Specifically, 
the AXP bursts we observed exhibit a wider range of durations, the correlation 
between burst fluence and duration is flatter than for SGRs, the observed AXP 
bursts are on average less energetic than observed SGR bursts, and the more 
energetic AXP bursts have the hardest spectra - the opposite of what is seen 
for SGRs. We conclude that the bursts are sufficiently similar that AXPs and 
SGRs can be considered united as a source class yet there are some interesting 
differences that may help determine what physically differentiates the two closely 
related manifestations of neutron stars. 
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1. INTRODUCTION 

Soft gamma repeaters (SGRs) are an exotic class of Galactic sources that are now 
commonly accepted as being magnetars - isolated, young neutron stars that are powered by 
the decay of an ultra-high magnetic field. The evidence for high surface fields (~ 10^^ — 
10^^ G) comes from several independent lines of reasoning (Duncan & Thompson 1992; 
Paczynski 1992; Thompson & Duncan 1995; Thompson & Duncan 1996). These include: 
the high dipolar magnetic fields implied by the spin properties of SGRs seen in quiescence 
under the assumption of magnetic dipole braking (Kouveliotou et al. 1998; Kouveliotou et al. 
1999); the requirement of a magnetar-strength field to confine the energy released in the tails 
of hyper-Eddington outbursts seen from two SGRs (Mazets et al. 1979; Hurley et al. 1999); 
the requirement of a high field to allow the decay rate necessary to power the burst and 
persistent emission (Thompson & Duncan 1996; Goldreich & Reisenegger 1992); and the 
magnetic suppression of the Thomson cross-section, which allows hyper-Eddington bursts to 
be observed (Paczynski 1992). For reviews of SGRs, see Kouvehotou (1999), Hurley (2000) 
and Thompson (2001). 

Anomalous X-ray pulsars (AXPs), another exotic class of Galactic neutron stars, have 
also been suggested to be magnetars (Thompson & Duncan 1996). This is because of their 
anomalously bright X-ray emission which can be explained neither by conventional binary 
accretion models nor rotation power (Mereghetti & Stella 1995). Also, their spin parameters, 
as for SGRs, imply large magnetic fields under standard assumptions of magnetic braking. 
They also have similar, though on average softer. X-ray spectra compared with those of 
SGRs in quiescence. However, unlike SGRs, in the > 20 yr since the discovery of the 
first AXP (Fahlman & Gregory 1981), none was seen to exhibit SGR-like bursts. For this 
reason, alternative models involving unconventional accretion scenarios have been proposed 
to explain AXP emission (van Paradijs et al. 1995; Chatterjee et al. 2000; Alpar 2001). See 
Israel et al. (2002) and Mereghetti et al. (2002) for reviews of AXPs. 

The magnetar model for AXPs was recently given a boost when SGR-like bursts were 
detected from two AXPs. Gavriil et al. (2002) reported on the discovery of two X-ray 
bursts in observations obtained in the direction of AXP IE 1048.1-5937. The temporal and 
spectral properties of those bursts were similar only to those seen only in SGRs. However, 
the AXP could not be definitely identified as the burster. On 2002 June 18, a major outburst 
was detected unambiguously from AXP IE 2259-1-586, involving over 80 bursts as well as 
significant spectral and timing changes in the persistent emission (Kaspi et al. 2003). Those 
bursts demonstrated that AXPs are capable of exhibiting behavior observed, until now, 
uniquely in SGRs, therefore implying a clear connection between the two source classes. 
Such a connection was predicted only by the magnetar model (Thompson & Duncan 1996). 
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However, the physical difference between the source classes is as yet unclear; Gavriil et al. 
(2002) and Kaspi et al. (2003) suggest that AXPs have higher surface magnetic fields than 
do SGRs, in spite of the evidence to the contrary from their spin-down properties. 

In this paper, we consider the statistical properties of the IE 2259+586 bursts in detail, 
in order to compare them quantitatively with SGR bursts, both to confirm that they have 
properties sufficiently similar that the two phenomena can definitely be unified, as well as to 
look for subtle differences that may offer clues regarding the physical distinction between the 
two classes. Statistical studies of magnetar bursts (e.g. Gogiis et al. 1999; Gogiis et al. 2000; 
Gogiis ct al. 2001) have the potential to yield important information regarding the burst 
energy injection and radiation mechanisms. Correlations between different burst properties, 
whether temporal and spectral, can be powerful model discriminators. Burst statistical 
properties can be compared with other physical phenomenon in order to assist in identifying 
their underlying cause; for example, they have been used to argue for important similarities 
between SGR bursts and earthquakes (Cheng et al. 1996). 

In this paper we present a comprehensive analysis of the properties of the bursts seen 
in the 2002 June 18 outburst of IE 2259+586. We present a study of the detailed outburst 
and post-outburst properties of the persistent and pulsed emission of IE 2259+586 in a 
companion paper (Woods et al. 2003). 

2. OBSERVATIONS AND ANALYSIS 

The results presented here were obtained using the Proportional Counter Array (PGA; 
Jahoda et al. 1996) on board the Rossi X-ray Timing Explorer (RXTE). The PGA consists of 
an array of five coUimated xenon/methane multi-anode proportional counter units (PCUs) 
operating in the 2-60 keV range, with a total effective area of approximately 6500 cm^ 
and a field of view of ~ 1° FWHM. We use RXTE to monitor all five known AXPs on 
a regular basis as part of a long-term monitoring campaign (see Gavriil & Kaspi 2002, 
and references therein). On 2002 June 18, during one of our regular monitoring observations 
{RXTE observation identification 70094-01-03-00) that commenced at UT 15:39:18, the AXP 
IE 2259+586 exhibited an SGR-hke outburst (see Fig. 1; Kaspi et al. 2003). The bursting 
behavior was detected by online RXTE monitors during the observation, and is clearly visible 
in the PGA "Standard 1" data. The observation spanned three orbits and had total on-source 
integration time 10.7 ks. Although some PCUs turned on/off during our observation, there 
were exactly three PCUs operational at all times. In addition to the standard data modes, 
data were collected in the GoodXenonwithPropane mode, which records the arrival time (with 
l-lis resolution) and energy (with 256-channel resolution) of every unrejected xenon event as 
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well as all the propane layer events. Processing of these data was done using software that 
operated directly on the the raw telemetry data. Photon arrival times were adjusted to the 
solar system barycenter using a source position of (J2000) RA 23'' OP" 08". 295, DEC +58° 
52' 44".45 (Patel et al. 2001) and the JPL DE200 planetary ephemeris. Note that following 
the outburst, Target of Opportunity observations of the source were initiated the next day 
and continued at different intervals over the subsequent weeks, however no more bursts were 
seen. 



To study the bursts quantitatively, we made use of the GoodXenonwithPropane data. 
Time scries were created separately for each PCU using all xenon layers. Light curves of 
various time bin widths (1/1024 s, 1/256 s, 1/64 s, 1/32 s and 1/16 s) were created to allow 
sensitivity to bursts on a range of time scales. The FTOOLs xtef ilt and meiketime were used 
to determine the intervals over which each PCU was off. We further restricted the data set 
by including only events in the energy range 2-20 keV. We used this energy range, which is 
larger than that used to study the quiescent pulsations (Gavriil Sz Kaspi 2002; Woods et al. 
2003), because of the much harder spectra of the bursts relative to the quiescent emission. 

The following procedure was performed separately for each PCU, in order to identify 
bursts. First, for each data set, the number of counts in the i^^ time bin was compared to 
a local mean //j. The local mean was calculated over a ~28 s (four pulse periods) stretch 
of data centered around the time bin being evaluated. A window of ~7 s (one pulse cycle) 
was also administered so that counts directly from, and immediately around, the point 
under investigation would not contribute to the local mean. During the outburst there was 
an increase in the pulsed flux (Kaspi et al. 2003; Woods et al. 2003), such that coherent 
pulsations were visible in our binned light curves. Because of this, for example, the apparent 
significance of bursts falling near a pulse peak would be artificially enhanced. To compensate 
for this effect, we first modelled the counts per time bin due to pulsations as: 



where A{(f), t) is the normalized amplitude of the pulsations as a function of pulse phase (p 
and time t. The parameters C and r are from an exponential fit to the pulsed flux evolution. 
We then calculated an adjusted local mean in the following way: 



j 

where the index j spans the windowed stretch of data used to calculate the local mean. 
For the number of counts in a time bin (rij) greater than the adjusted local mean (Aj), the 



2.1. The Burst Identification Algorithm 



(1) 




(2) 
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probability of those counts occurring by random chance is given by 



n,;! 



(3) 



As the probabihty for each PCU is independent, we calculated the total probability {Pm) 
of observing a burst simultaneously by all operational PCUs as 



Pi,toi^l[Pi,k, (4) 



k=0 

and k corresponds to the PCU under consideration. If a particular PCU were inoperable we 
set Pi^k — 1- Events which registered a value of Pj^tot < 0.01/ N, where N is the total number 
of time bins searched, were flagged as bursts, and were subject to further investigation. 

The significance of the number of counts in a time bin can be underestimated if there 
are one or more bursts in the interval used as the local mean. For this reason, once a burst 
was identified it was removed from the light curve, and the burst identifying procedure was 
repeated until there were no additional bursts returned. 



3. RESULTS 

3.1. Burst Statistics 

Our burst searching algorithm returned 80 significant bursts from the 2002 June 18 
observation. The number of bursts identified depended on the time resolution used: 26%, 
55%, 76%, 83% and 74% of all identified bursts were fiaggcd at 1/1024 s, 1/256 s, 1/64 s, 
1/32 s and 1/16 s time resolution, respectively. Most bursts were single-peaked and had 
durations <1 s. A small handful (~12) were bright and had clear fast- rise, exponential decay 
morphology. Four bursts were multi-peaked. A variety of burst morphologies is shown in 
Figure 2. Some bursts (~5%) were approximately symmetric, a few (~3%) fell faster than 
they rose while most fell slower than they rose (see §3.1.4). 



3.1.1. Burst Event Times and Phase 

The time of each burst was initially defined, using binned light curves, to be the midpoint 
of the bin having the most counts. To increase the precision of the burst time we refined this 
value, using the event data which comprised this time bin, to be the midpoint of the times 
of the events having the smallest temporal separation. We also calculated the occurance in 



-6- 



pulse phase for each burst using the time of the burst peak and the rotational ephemeris 
given by Kaspi et al. (2003). Comparing the burst phase distribution to the pulse profile of 
IE 2259+586 at the time of the outburst, a correlation is seen (Fig. 3), where most of the 
bursts tend to occur when the pulsed intensity is high. We note that the two bursts seen 
from the AXP IE 1048.1-5937 (Gavriil et al. 2002) were also coincident with the pulse peak, 
which strengthens the argument that IE 1048.1-5937 was the source of those bursts. 



3.1.2. Burst Durations and Fluence 

The Tgo duration is the time between when 5% and 95% of the total background- 
subtracted burst counts have been accumulated (e.g. Gogiis et al. 2001). The background 
count rate was determined by averaging a hand-selected burst-free region before and after the 
burst. This typically consisted of two intervals of 1 s before and after the burst in question. 
The integrated background-subtracted counts were then fit to a step function plus a linear 
term using least-squares fitting. The height of the step-function corresponds to the total 
burst fluence F (in counts) and the slope of the line corresponds to any background counts 
that were improperly subtracted. 

SGR Tgo distributions follow a log-normal distribution, defined as 

P(T9o, /X, <t) = ^--= exp 

log(Tv27r 

whose mean and standard deviation vary with source (e.g. Gogii§ et al. 2001). At first we 
fit the measured values of Tgo for the IE 2259+586 bursts with this model and found it to 
characterize the distribution well. In Equation 5 the parameters log/x and log a correspond 
to the mean and standard deviation of the the log Tgo values. The mean of the Tgo values 
is given by n and the range for one standard deviation corresponds to {fj,a~^ , fj^a) . The 
best-fit II and a were determined by maximum likelihood testing. The latter allowed us 
to extract model parameters that are independent of the arbitrarily chosen histogram bin 
widths. Specifically, the best-fit parameters were those which maximize the statistic 

N 

M = J2^ogP{T,o,i,l^,a), (6) 

i=l 

where is the number of bursts. Figure 4 shows the distribution, and best-fit log-normal 
model for the measured values. We found that our Tgo distribution has mean fi — 97.9 ms 
with a range of 18.2-527.2 ms for one standard deviation. Note however that for low signal- 
to-noise bursts, Tgo can be substantially underestimated. We describe how we corrected for 
this problem and obtained slightly modified best-fit log-normal parameters in §3.1.5 below. 



1 / log Tgo - \0gfl 

2 V loga 



(5) 
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The fluences measured as described above were then grouped in equispaced logarithmic 
bins. The distribution of burst fluences is displayed in Figure 5. The low-end fluences are 
underrepresented because of sensitivity drop-off. Excluding the points having fluence <20 
PCA counts, the distribution is well modeled by a simple power law. Using least-squares 
fitting we find a best-fit power-law index of —0.7 ± 0.1, which corresponds to a differential 
spectrum dN/dF oc ir-i-7±o.i pj-om the plot, it is clear that the fluences span approximately 
two orders of magnitude. For our calibration of the fluences in CGS units, sec §3.2.4. 

Gogii§ et al. (2001) also find a clear correlation between burst durations and total burst 
fiucncc. In Figure 6, wc plot fiucncc versus T90. A correlation can clearly be seen. To 
quantify it, we grouped the T90 values in equispaced logarithmic bins and determined group- 
averaged fluences for each bin. Least-squares fltting to a simple power-law model yields 



Burst peak fluxes were determined from the event data using the following algorithm. A 
box-car integrator of width 62.5 ms was translated through the event data. The procedure 
began and ended when the center of the box-car was at half a box-car width before and 
after the time of the burst peak (as determined in § 3.1.1). At each box-car step a flux 
measurement was made by integrating the number of events and dividing by the box-car 
width. The burst peak flux was assigned the largest such flux measurement. We then 
grouped our peak fluxes in equispaced logarithmic bins. The distribution of peak fluxes is 
shown in Figure 7. 

Our burst-identifying algorithm is less sensitive to bursts of smaller peak flux. To 
compensate for this effect, we ran the following simulation. We took a hand-selected 1- 
ks long burst-free region from our observed IE 2259+586 light curve binned with 62.5-ms 
resolution. We then injected a simulated burst having peak flux fp at a random position in 
the light curve. We modelled the burst by a top-hat function of width 62.5 ms (one time bin) 
and height fp x 62.5 ms. We then ran our burst-identifying algorithm as described in § 2.1. 
We repeated this procedure for Ni iterations and determined Ng, the number of successful 
burst identiflcations for that simulated peak flux. We repeated the procedure for various 
peak fluxes and determined the probabihty of detecting a burst P — Ng/ as a function of 
peak flux fp. We found that P could be well modelled by the following analytic function 



with /o = 309.84 cts s ^ and k — 58.21 cts s ^. We then used this function to correct 



F oc T+°•5^±°•°^ with reduced = 1-0. 



3.1.3. Burst Peak Fluxes 




(7) 
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our peak flux distribution (see Fig. 7, boxes) . Using least-squares fltting we found that the 
corrected distribution is weU modelled by a simple power law with index —1.42 ± 0.13. For 
our calibration of these peak fluxes in CGS units, see §3.2.4. 



Burst rise and fall times were obtained from the event data by maximizing the likelihood 
of the assumed probability distribution 



where B represents the background count rate, Cp represents the background-subtracted 
count rate at the time of the burst peak tp, and tr and tf represent the burst rise and fall 
times, respectively. The parameter A is a normalizing factor ensuring unit probability over 
the interval of interest. This model characterized the bursts well - see the left panels of 
Figure 2 (dotted line) for examples. Burst rise and fall time distributions are displayed in 
Figure 8, with best-flt log-normal models determined via maximum-likelihood testing. For 
the rise time distribution, we flnd a mean of 2.43 ms and a range of 0.51-11.51 ms for one 
standard deviation, with reduced — For the fall time distribution, we flnd mean 
13.21 ms and a range of 3.52-49.55 ms for one standard deviation, and a reduced = 0-2. 
In order to better quantify burst morphologies we also show the ratio of burst rise times to 
fall times [tr/tf] Fig. 8). On average, bursts rise faster than they fall, however this is not 
universally true. Again fitting a log-normal distribution, wc find mean 0.18 and a range of 
0.03-1.08 for one standard deviation, with reduced = 3.7. The latter fit is poor because 
the distribution is clearly skewed toward shorter rise times. The asymmetry of the typical 
burst can also be seen in Figure 9, where the distribution of tr/Tgo is plotted. 



G5gii§ et al. (2001) showed that in the low signal-to-noise regime, the value of Tgo can 
be underestimated. To account for this, a model light curve was generated for each burst, 
having the form of Equation 8. Peak flux, rise time and fall time were flxed at the values 
measured for that particular burst. The simulated light curve was then integrated and the 
model duration (Tgo,^) was measured by the same procedure outlined in §3.1.2. We then 
repeated the procedure with noise added to the simulated light curve. The noise was drawn 
from a Poissonian distribution having mean equal to the measured background rate of the 



3.1.4- Burst Rise Times and Fall Times 




(8) 



3.1.5. Corrected Tgo Values 
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burst under investigation. We repeated the procedure for 200 realizations of noise. For each 
iteration (i) wc measured the duration (Tgo^i). The simulated durations (Tgo^i) were normally 
distributed and the mean of this distribution (Tgo^^ ) allowed us to calculate a correction 
factor J-'V = 1 — Tgo^m/Tgo^s. The corrected Tgo distribution is shown in Figure 4. The 
best-fit mean is 99.31 ms with a range of 14.4-683.9 ms for one standard deviation. 

3.1.6. Burst Waiting Times 

SGR waiting times (AT), defined as the temporal separations of adjacent bursts, are 
found to follow log- normal distributions (Gogiis et al. 1999; Gogiis et al. 2000). We mea- 
sured the waiting time for the IE 2259+586 events, excluding those interrupted by Earth 
occultations. Figure 10 displays our AT distribution with the best-fit log-normal model as 
determined by maximum likelihood testing. The best-fit parameters are mean of 46.7 s and 
a range of 10.5-208.4 s for one standard deviation, with reduced ~ 0.6. We find no 
correlation between the burst energy, duration and the waiting time until the next burst, 
nor with the elapsed time since the previous burst. 

Note however that the burst rate clearly decreased during the observation (see Fig 1). 
This is made clear by the bottom panel of Figure 10 which shows a correlation between the 
waiting time (AT) and the burst peak time (tp). We fit this correlation to a power-law model 
using least-squares fitting, which reveals that AT = 0.11 x tp^-^^ . This correlation implies 
that the mean of our waiting time distribution depends on the time at which we started 
observing the outburst. 

3.2. Burst Spectroscopy 

3.2.1. Individual Burst Spectra 

Spectra for each burst were extracted with the 256 spectral bins over the PGA range 
grouped by a factor of 4 in order to increase the signal-to-noise ratio per spectral bin. The 
same background intervals selected in measuring Tgg were used in the spectral analysis (see 
§3.1.2). In all spectral analyses, energies below 2 keV and above 60 keV were ignored, 
leaving on average 33 spectral channels for fitting. The regrouped spectra along with their 
background estimators were used as input to the X-ray spectral fitting software package 
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XSPEC^. Response matrices were created using the FTOOLs xtef ilt and pcarsp. We fit the 

28 most fluent bursts with a photoelectrically absorbed power law of index F, holding only 
Nh flxed at 0.93 x 10^^ cm~^ (the value found by Patel et al. 2001). The distribution of 
spectral indices is shown in Figure 11. We find a mean spectral index of F = 1.35 with 
standard deviation 0.43. 



3.2.2. Hardness Ratios 

Gogiis et al. (2001) noted that SGR bursts tend to soften with increasing burst energy. 
We studied the hardness ratio/fluence relationship by extracting spectra and creating re- 
sponse matrices separately for each burst. Hardness ratios were defined as the ratio of the 
counts in the 10-60 keV band to those in the 2-10 kcV band as in Gogiis et al. (2001). 
Also following Gogii§ et al. (2001), we divided the bursts into equispaced logarithmic fiu- 
ence bins and calculated a weighted average hardness ratio for each bin. Figure 12 shows 
the weighted mean hardness ratios as a function of fluence. A clear positive correlation is 
seen. We repeated the procedure for different definitions of hardness ratio and found similar 
correlations. We further confirmed this trend by considering the 28 most fiuent bursts for 
which spectral indexes F could be reliably and precisely constrained. All had F well below 
the mean value. 



3.2.3. Absence of Spectral Lines and the Average Burst Spectrum 

Possible spectral features have been reported in a burst from the AXP IE 1048.1- 
5937 (Gavriil et al. 2002) and from bursts from two SGRs (Strohmayer & Ibrahim 2000; 
Ibrahim et al. 2002, 2003). In no spectrum of any burst for IE 2259+586 did we detect a 
significant feature. In order to amplify any low-level spectral feature common to all bursts, 
we combined individual burst spectra to create a grand average spectrum. We summed the 
burst and background spectra described in the previous section using the FTOOL sumpha. 
Response matrices were scaled and added using the FTOOL addpha. In order to search for 
features in the residuals, we fit the combined spectrum to a simple photoelectrically absorbed 
power law. The residuals showed no evidence of significant spectral features. 



^http: / /xspec. gsfc.nasa.gov 
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3.2.4- Calibrating Fluence and Flux 

Determining peak flux and total fluence distributions in CGS units requires spectral 
fitting. However most bursts were too faint to allow spectral parameters to be determined 
with interesting precision. The problem was worse for the peak fluxes since even the brighter 
bursts generally had too few counts to meaningfully constrain the spectrum. Therefore, we 
devised an alternate way of converting between PCA counts and CGS units. We took the 
spectra of the 40 most luminous bursts extracted over their Tgo duration and fit them with 
photoelectrically absorbed power laws. However this time, for consistency, we held F fixed at 
the mean of our spectral index distribution. We multiplied the flux (in units of erg s~^ cm~^) 
in the 2-60 keV range returned by the fit by its respective Tgo duration to obtain a fluence 
in erg cm~^. We then considered the 2-60 keV fluence in counts as determined in §3.1.2 as a 
function of the fluence in CGS units and determined the proportionality constant between the 
two using least-squares fitting. This constant was found to be 8.226 x 10~^^ erg cm~^ cts~^. 
In §3.2.2 we found significant spectral evolution as a function of fluence. A change of la in 
spectral index F corresponds to a change by a factor of ~ 1.5 in our cahbration constant. The 
same procedure and constant apphes for the peak fluxes. The CGS energy scales are shown 
at the top of Figures 5 and 7. The fluences in the 2~60 keV band range from ~ 5 x 10~^^ to 
~ 7 X 10^^ erg cm~^. These imply burst energies in the range ~ 5 x 10'^'^ to ~ 7 x 10'^^ erg, 
assuming isotropic emission and a distance of 3 kpc to the source (Kothes et al. 2002). The 
sum total of all burst fluences is 5.6 x 10^^ erg cm~^, corresponding to energy 6.0 x 10'^'' erg (2- 
60 keV). Peak fluxes in a 61.25-ms time bin range from ~ 1 x 10~^ to ~ 1 x 10~^ erg cm~^ s^^, 
which imply peak luminosites in the range ~ 1 x 10^^ to ~ 1 x 10^^ erg s~^. On shorter time 
scales we find 5 bursts with peak fluxes which are super-Eddington. The peak fluxes in a 
1/2048 s time bin for these bursts range from ~ 2 x 10^^ to ~ 8 x 10^^ erg s~^. 

4. DISCUSSION 

As we describe below, many of the properties of the bursts seen from IE 2259+586 during 
its 2002 June 18 outburst are very similar to those seen in SGRs 1806—20 and 1900+14. 
However, there are some quantitative differences. Next we compare the various measured 
quantities for the AXP and SGR bursts. Note that our comparisons focus primarily on PCA 
observations of SGRs for consistency of spectral and temporal response. 

The mean Tgo value of 99.31 ms (see §3.1.5 and Fig. 4) is very similar to those seen 
for SGRs 1806-20 and 1900+14: 161.8 ms and 93.9 ms, respectively. Gogii§ et al. (2001) 
suggested that the difference between these values for the two SGRs is a result of a different 
intrinsic physical property of the sources, such as the strength of magnetic fleld, or the size of 
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the active region. Given the generally softer persistent emission spectra of AXPs compared 
to SGRs, as well as the less frequent outbursts of the AXPs, it is reasonable to suspect that 
the two source classes differ also by some physical property; age (Kouveliotou et al. 1998; 
Gaensler et al. 2001), magnetic field (Gavriil et al. 2002; Kaspi et al. 2003) and progenitor 
mass (Gaensler et al. 2001) have been proposed. The similarity of the burst durations of all 
three sources implies, however, that the physical property resulting in different mean burst 
durations must be different from that which results in different average spectra and outburst 
frequency. 

The standard deviation of the Tgo distribution for the IE 2259+586 bursts is much 
larger than is seen for the SGR bursts. For IE 2259+586, the la range is from ~14 ms to 
~684 ms or 1.7 magnitudes. For SGRs 1806—20 and 1900+14, the corresponding range in 
durations is 0.68 and 0.70 magnitudes. The lower bound on the IE 2259+586 distribution 
may be artificially lower due to the shorter time scales searched in this work as compared to 
Gogii§ et al. (2001) who searched for SGR bursts on the 0.125 s time scale. However, such a 
wide range of durations is seen even when faint bursts are omitted from the Tgo distribution 
of IE 2259+586. Gogiis et al. (2001) argued that if the "trapped fireball" model, which 
describes the giant SGR bursts well, also applies to the fainter bursts, then the narrowness 
of the Tgo distribution compared with the wide range of fiuences demands a planar fireball 
geometry. This is because the duration of the burst is hmited by the rate of coohng through 
the radiative fireball surface layer. For IE 2259+586, the Tgo range is larger than the 
fluence range, indicating that if the fireball model applies, a planar fireball geometry is not 
supported. 

The distribution of burst fiuences for IE 2259+586 is remarkably similar to those 
seen in SGRs. For the IE 2259+586 bursts, we find a fiuence distribution dN/dF oc 
^-i.7±o.i (pjg 5) Gogiis et al. (2000) showed that for the PGA, the fiuence distribution for 
SGR 1806—20 is well described by a power law of index —1.43 ± 0.06, while at higher burst 
energies, the index steepens to —1.7. For SGR 1900+14, Gogiis et al. (1999) found an index 
of — 1.66lo}2 extending over the full range of burst fiuences. The good agreement of the 
fluence distribution indices shows that for a given outburst intensity (i.e. the normalization 
of the fluence distribution), the average burst energy is the same for IE 2259+586 as it is 
for these two SGRs. The difference between the SGR outbursts that are routinely detected 
by IPN detectors and this outburst from IE 2259+586 which was not detected by the IPN 
is the SGR outbursts have shown higher outburst intensities. Since we know that the SGRs 
spend most of their time in quiescence when the fiuence distribution normalization is zero 
(or near zero), the dynamic range of the outburst intensities in SGRs is larger than has been 
observed thus far in IE 2259+586. This difference in range is intrinsically even larger when 
one considers that IE 2259+586 is believed to be significantly closer (3 kpc) than either of 
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these two SGRs (~15 kpc, Vrba et al. 2000; Corbel et al. 1997). 

Cheng et al. (1996) noted the similarity of the fluence distribution index for SGR 1806—20 
with that determined empirically for earthquakes (Gutenberg & Richter 1956a,b, 1965), and 
also for the distribution of earthquake energies found in computer simulations (Katz 1986). 
However, solar flares also show a size distribution with exponents ranging from 1.53 to 1.73 
(Crosby et al. 1993; Lu et al. 1993). Magnetars are not clearly physically analogous to either 
system; in magnetars, magnetic stresses are thought to result in stellar crust cracking, which 
is not the case for earthquakes. The bursts could be magnetic reconnections as in solar 
flares (Lyutikov 2002), however in the solar case there is no solid crust to yield, unlike in 
magnetars. The similarity of the distributions could be explained as being a result of the 
phenomena of self-organized criticality (Bak et al. 1988), in which a system is dynamically 
attracted (i.e. self-organized) to a critical, spatially self-similar state which is just barely 
stable to perturbations. In other words, the burst statistics alone do not constrain their 
physical origin. 

It is not possible to compare peak flux distributions as none are published for SGRs. 
For the AXP, the range of 2-60 keV peak flux for the 62.5-ms time scale spans a factor of 
~100, ranging from ~ 1 x 10~^ to ~ 1 x 10~^ erg cm~^ s~^. which, for a distance of 3 kpc, 

corresponds to luminosities of ~ 1 x 10^^ to ~ 1 x 10"^^ erg s~^. At time scales as short as 
1/2048 s we find peak fiuxes as high as ~ 8 x 10^® erg s~^. Thus 5 bursts are above the 
Eddington limit on this time scale. 

As in SGRs, the fluences of the IE 2259-1-586 bursts are significantly positively correlated 
with Tgo (Fig. 6). However there is one difference: for the AXP, the relationship is well 
described by a power law of index +0.54 ± 0.08, while for SGRs 1806—20 and 1900+14, 
Gogii§ et al. (2001) found +1.05 ± 0.16 and +0.91 ± 0.07, respectively. Thus the power-law 
index for AXPs is half that seen in SGRs. It is important to recognize, however, that severe 
selection effects are at work here. Specifically, as discussed in §3.1.2, we are less sensitive 
to low-fluence bursts. This is particularly true for bursts having long rise times, which will 
tend to have long Tgo values. Thus there are severe selection effects against finding bursts 
in the bottom right-hand portion of Figure 6, as there are in similar analyses for SGRs. 
Therefore the above correlation should really be seen as an upper envelope to the phase 
space available to the burst. By contrast, our sensitivity to bursts that would sit in the 
upper left-hand corner of the plot is generally enhanced relative to the populated region, 
indicating the absence of bursts in this part of phase space is genuine. 

The morphologies of the AXP and SGR bursts are similar, with most being asymmetric, 
with faster rises than decays. Rise and fall time distributions for the SGRs have not been 
published, so we cannot compare those parameters directly, nor the ratio of the two. Gogii§ 
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et al. (2001) showed the distribution of the ratio i^/Tgo for SGRs 1806-20 and 1900+14; 
the same plot for IE 2259+586 fooks similar (Fig. 9). 

The waiting time distributions of the AXP and SGRs are very similar. All are well 
described by log-normal distributions. This is similar to what is seen in other self-organized 
critical systems, such as earthquakes (Nishenko & Buland 1987). For IE 2259+586, we find 

a mean waiting time between bursts of 47 s, and range of 10-208 s. Gogiis et al. (1999) found 
~49 s for SGR 1900+14, and G6giis et al. (2000) found ~97 s for SGR 1806-20, with range 
between ~0.1 and 1000 s for both, very similar to our results. The absence of correlation of 
waiting time and burst fiuence for the AXP is similar to that seen for SGRs (Gogiis et al. 
1999; Gogiis et al. 2000), although Gogiis et al. (1999) report an ant icor relation between 
time since the previous burst and burst energy. We do not see this for the AXP, nor do 
Gogiis et al. (2000) observe it for SGR 1806-20. 

One striking difference between the AXP and SGR bursts is in the relationship between 
spectral hardness ratio and fiuence. For SGR 1806—20, Gogii§ et al. (2001) found that the 
more energetic bursts are spectrally softer, regardless of burst morphology. This was not 
seen for SGR 1900+14, however. Our analysis (see Fig. 12) shows the opposite behavior to 
that seen in SGR 1806—20, with the more energetic bursts having harder spectra. Gogii§ 
et al. (2001) argued that the behavior seen for SGRs could be explained either by the 
emitting plasma being in local thermodynamic equilibrium, having radiative area decreasing 
for lower fiuences, or by the spectral intensity of the radiation field being below that of a 
blackbody, hence the emitting plasma temperature T remaining in a narrow range, being 
higher at lower luminosities. Which of these two applies depends on the rate of energy 
injection into the magnetosphere; the latter applies only if the luminosity is less than ~ 
10^^(y^/^/10 km) erg s~^ where V is the injection region, assuming a spherical geometry. 
Glearly neither can apply for the AXP. Gogii§ et al. (2001) imply that blackbody emission 
from a constant radius predicts the relationship between hardness and fiuence that we find 
for the AXP. However for the AXP, naively taking Figure 6 at face value, F oc Tg5^. Hence 
La OC F~^, so blackbody emission from a constant radius predicts T oc the opposite 

to what we have observed. We note further that the range of hardness ratios for the AXP 
bursts is slightly greater than it is for the SGRs. For IE 2259+586, hardness ratios (for 
bursts having 10^-10^ counts) range from ~0. 54-0. 85, while the range is ~0. 82-0. 95 for 
SGR 1806-20, and ~ 0.63-0.67 for SGR 1900+14 (Gogii§ et al. 2001). It should be noted 
however that we identified bursts (see §2.1) using a different energy range (2-20 keV) than 
Gogii§ et al. (2001), who used the full bandpass of the PGA. This would make us more 
sensitive to softer bursts which would affect the dynamic range of the hardness ratios we 
measured. Perhaps interestingly, for the SGRs, F oc Tgo, so the La = F/Tgo ~ constant, 
and for constant radiative area and blackbody emission, one expects T ~ constant, closer 
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to what is observed for SGRs than for AXPs. Thus, although blackbody emission from 

a constant radius (not surprisingly) does not describe any of the data well, it does seem 
possible that the flatter dependence of fluence on Tgg, the inverted dependence of hardness 
on fluence relative to the SGRs, and the greater range of hardness in the AXP bursts may all 
be related phenomena telling us something interesting about the physical distinction between 
these closely related sources. 

We have stated that outbursts from AXPs similar to or larger than the one studied 
here are less frequent than are those from SGRs. Of course, given that we have observed 
only one AXP outburst, and that this outburst was energetically smaller and fainter than 
observed SGR outbursts, making a meaningful comparison of their outburst rate is very 
difficult. We can estimate the rate of AXP outbursts of the magnitude of the 2002 June 
18 event as follows. We consider data from only our RXTE PGA monitoring program, 
as it provides a consistent quasi-regular ly sampled data set with a single instrument. The 
monitoring program for IE 2259+586 has extended over nearly 7 yr with only one such 
outburst detected; even though the bursting appears to have been relatively short-lived, the 
effects of a glitch of even much smaller size would easily have been detected throughout 
the data span. We make the admittedly speculative assumption that all such outbursts 
are accompanied by comparably sized glitches. A comparable glitch in AXP IRXS 1708- 
4009 was recently detected in 5.4 yr of monitoring without evidence for radiative outburst, 
however the sparse observations could have missed one (Kaspi & Gavriil 2003; DaU'Gsso 
et al. 2003). Two small bursts have been seen in 6.8 yr of timing of AXP IE 1048.1-5937 
(Gavriil et al. 2002), and its timing behavior suggests that many ghtches could be occurring 
(Kaspi et al. 2001), however no other evidence for radiative outbursts has been found. No 
activity of any kind, apart from apparently simple timing noise, has been seen in 6.5 yr of 
timing of 4U 0142+61 (Gavriil & Kaspi 2002) or in 4.3 yr of timing IE 1841-045 (Gotthelf 
et al. 2002). If we omit IE 1048.1-5937 whose timing behavior we do not fully understand, 
we can estimate a rough AXP outburst rate of one every 11 yr, assuming that the glitch 
in IRXS 1708-4009 was indeed a similar outburst, or one every ~22 yr if not. SGRs, by 
contrast, burst much more frequently, reach higher intensities, and persist for longer periods 
of time. The monitoring of the SGRs with the RXTE PGA has not been as regular as 
for the AXPs due to less optimal observing conditions for the SGRs (lower pulsed fractions, 
source flux, stronger timing noise, etc.), therefore, we cannot make a direct comparison of the 
outburst recurrence rate using the PGA data. We can, however, make a rough estimate of 
the recurrence rate using results obtained with the Burst and Transient Source Experiment 
(BATSE) that flew aboard the Compton Gamma-Ray Observatory. The advantage of using 
BATSE to estimate the SGR outburst rate is its uniform and dense coverage in time due to 
its "all-sky" FOV. The disadvantage is that BATSE is much less sensitive to SGR bursts than 
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is the PCA (e.g. Gogiis et al. 1999). Since SGR/AXP burst energies follow a steep power- 
law distribution, the outburst recurrence rate is a strong function of detector sensitivity. It 
follows that an outburst recurrence rate determined by BATSE will then be a lower limit 
to the rate for the more sensitive PCA. Moreover, the relative distances of AXPs and SGRs 
must be considered when determining intrinsic source rates for a given luminosity or total 
energy as opposed to peak flux and fluence. With these factors in mind, we now estimate the 
SGR outburst recurrence rate at the BATSE sensitivity level. BATSE was in operation for 
9.1 yr from 1991 April through 2000 June. During that time, three of the four known SGRs 
entered outburst (Kouvehotou et al. 1993; Kouveliotou et al. 1994; Woods et al. 1999; Gogii§ 
et al. 2001), some multiple times. Here, we define an outburst as a collection of bursts (i.e. 
more than two) where the separation between consecutive bursts never exceeds one month. 
Using the results reported in Gogii§ et al. (2001), the number of SGR outbursts detected 
during this time interval is 14. This yields an outburst rate for the SGRs of once every ~2.6 
years. Recall, this is a lower limit to the rate at the PCA sensitivity level. Thus the SGRs 
clearly undergo outbursts more frequently than do AXPs. 



5. CONCLUSIONS 

The bursts we have observed for IE 2259+586 are clearly similar to those seen uniquely 
in SGRs. As concluded by Gavriil et al. (2002) and Kaspi et al. (2003), AXPs and SGRs 
clearly share a common nature, as has been predicted by the magnetar model. In this 
paper, we have done a quantitative analysis of the IE 2259+586 bursts seen on 2002 June 
18, and compared our results with those obtained for the two best-studied SGRs, 1806—20 
and 1900+14. We summarize our results as follows. The bursts seen in the 2002 June 18 
outburst of IE 2259+586 are qualitatively similar to those seen in SGRs, and in many ways 
quantitatively similar. Specifically: 

• the mean burst durations are similar 

• the differential burst fiuence spectrum is well described by a power law of index —1.7, 
similar to those seen in SGRs (and earthquakes and solar fiares) 

• burst fluences are positively correlated with burst durations 

• the distribution of and mean waiting times are similar 

• the burst morphologies are generally asymmetric, with rise times usually shorter than 
burst durations 
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However, there are some interesting quantitative differences between the properties of the 
AXP and SGR bursts. These may help shed hght on the physical difference(s) between these 
classes. The differences can be summarized as: 

• there is a significant correlation of burst phase with pulsed intensity, unlike in SGRs 

• the AXP bursts have a wider range of burst duration (though this may be partly due 
to different analyses procedures) 

• the correlation of burst fluence with duration is flatter for AXPs than it is for SGRs 
(although when selection effects are considered, this correlation should really be seen 
as an upper envelope for AXPs and SGRs) 

• the fluences for the AXP bursts are generally smaller than are in observed SGR bursts 

• the more energetic AXP bursts have the hardest spectra, whereas for SGR bursts, they 
have the softest spectra 

• under reasonable assumptions, SGRs undergo outbursts much more frequently than 
do AXPs 

Given the rarity of AXP bursts coupled with the unique information that detection of 
such bursts provides, observing more outbursts is obviously desirable. Continued monitoring 
is thus clearly warranted, and RXTE with its large area and flexible scheduling is the obvious 
instrument of choice. 

We are grateful to C. Kouveliotou, M. Lyutikov, S. Ransom, M.S.E. Roberts, D. Smith, 
and C. Thompson for useful discussions. This work was supported in part by NSERC, 
NATEQ, CIAR and NASA. This research has made use of data obtained through the 
High Energy Astrophysics Science Archive Research Center Online Service, provided by 
the NASA/Goddard Space Flight Center. 
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Fig. 1.— 2-60 keV RXTE /PC A light curve for IE 2259+586 on 2002 June 18, at 62.5-ms 
resolution. The gaps are Earth occultations. 
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Fig. 2. — Three different examples of bursts seen in the 2002 June 18 outburst of 
IE 2259+586. Left: Sample background-subtracted hght curves in the energy range 2- 
60 keV with 1/32 s (top), 1/512 s (middle) and 1/2048 s (bottom) time resolution. The 
dotted line shows the model fit to the data in order to measure burst rise and fall times (see 
§3.1.4 for details). Right: Cumulative background-subtracted counts for each burst. The 
vertical dotted line shows the location of the burst peak. The horizontal dotted line shows 
the level used in determining the burst fluence. See §3.1.2 for details. 
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Fig. 3. — Distribution of the pulse phases of IE 2259+586 which correspond to the times of 
the burst peaks (solid points). The solid curve is the folded 2-60 keV light curve of the 2002 
June 18 observation with the bursts omitted. 
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Fig. 4. — Distribution of Tgo durations for the bursts observed from IE 2259+586. The 
sohd histogram hne shows the observed binned distribution (see §3.1.2), while the dashed 
histogram hne shows the corrected distribution (see §3.1.5). The solid curve represents 
the best-fit log-normal model for the observed data, as determined by maximum-likehhood 
testing. The dashed curve is the best-fit log-normal model for the corrected data. This fit 
has mean 99.31 ms and standard deviation of a factor of 6.9. 
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Fig. 5. — Distribution of the 2-60 keV fluence F for each burst observed from IE 2259+586. 
Sohd points represent average values of fluence in equispaced logarithmic bins for which our 
observations had full sensitivity. The open points suffered from reduced sensitivity. The 
best-fit line was determined using the solid points only and is shown as a solid line; the 
dashed fines are its extrapolation. The slope of this line is — 0.7 ± 0.1, which corresponds to 
dN/dF oc F-i-^. 



-26- 




Fig. 6. — Burst 2-60 keV fluence versus Tgo. The open points represent individual bursts. 
The sohd points represent binned averages. The sohd hne represents the best-fit power law 
for the binned averages. The slope of the line is +0.54 ± 0.08. 
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Fig. 7. — Distribution of burst peak flux for 62.5-ms time binning. The diamonds are 
observed averages in equispaced logarithmic bins. Our sensitivity is significantly reduced at 
low peak fluxes. The corrected values, determined using simulations described in § 3.1.3 are 
shown by open squares. The corrected flux bins were fit with a power law, shown by a line. 
The slope is -1.42 ±0.13. 
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Fig. 8. — Distribution of burst rise (t^) and fall (t/) times (see §3.1.4). Bottom left: Dis- 
tribution of fall times t^- Bottom right: Distribution of fall times tf. Top: Distribution of 
tr/tf. In all cases, the solid line represents the best fit log- normal model, as determined by 
maximum-likelihood testing. 
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Fig. 9. — Distribution of the ratio of burst rise time tr to duration Tqq. 
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Fig. 10. — Top: Distribution of the waiting time between successive bursts. The sohd Une 
represents the best fit log-normal model, as determined by maximum-hkehhood testing. The 
mean is 46.8 s, and standard deviation of a factor 4.4. Bottom: Waiting time as a function 
of event time. The line represents the best-fit power law model. 
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Fig. 11. — Distribution of spectral indices (F) for the 28 most fluent bursts. See §3.2.1 
for details. The curve is the best-fit gaussian model. This fit has mean 1.35 and standard 
deviation 0.43. 
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Fig. 12. — Hardness ratio (if) versus fluence (F). Hardness ratio is deflned as the ratio of the 
number of PCA counts in the 10-60 keV band to that in the 2-10 keV band. The open points 
are hardness ratio measurements for individual bursts. The sohd points are weighted averages 
of hardness ratios for bursts in equispaced logarithmic fluence bins. The line represents the 
best-fit logarithmic function for the weighted averages, H = 0.31 x logF — 0.09. 



